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We study the production of the hghtest neutraUnos in the radiative process e"'"e~ — > XiXi7 in 
supersymmetric models with grand unification. We consider models wherein the standard model 
gauge group SU{'i)c x SU{2)l x U{1)y is unified into the grand unified gauge groups SU{5), or 
50(10) . We study this process at energies that may be accessible at a future International Linear 
Collider. We compare and contrast the dependence of the signal cross section on the grand unified 
gauge group, and different representations of the grand unified gauge group, into which the standard 
model gauge group is unified. We carry out a comprehensive study of the radiative production 
process which includes higher order QED corrections in our calculations. In addition we carry out 
a detailed study of the background to the signal process coming from the Standard Model radiative 
^SJ , neutrino production e"'"e~ — >■ as well as from the radiative production of the scalar partners of 

the neutrinos (sneutrinos) e^e~ — > i>i>*j. The latter can be a major supersymmetric background 
' to the radiative production of neutralinos when the sneutrinos decay invisibly. It is likely that the 

, radiative production of the lightest neutralinos may be a viable channel to study supersymmetric 

partners of the Standard Model particles at the first stage of a International Linear Collider, where 
heavier sparticles may be too heavy to be produced in pairs. 
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PACS numbers: ll.30.Pb, 12.60.Jv, 14.80.Ly 



I. INTRODUCTION 

^: 

Supersymmetry (SUSY) [l| is a leading candidate for physics beyond the standard model (SM). In super- 
symmetric models the Higgs sector of the SM is technically natural ■ Since supersymmetry is not an exact 
^ • symmetry in nature, it must be broken in realistic models of supersymmetry. Although the precise manner in 
which SUSY is broken is not known at present, in actual practice necessary SUSY breaking can be introduced 
through soft supersymmetry breaking terms that do not reintroduce quadratic divergences in the Higgs mass, 
and thereby do not disturb the stability of the hierarchy between the weak scale and a large scale (grand 
^ , unified (GUT), or Planck scale). The simplest implementation of the idea of softly broken supersymmetry is 
the Minimal Supersymmetric Standard Model (MSSM) obtained by introducing the supersymmetric partners 
of the SM states, and introducing an additional Higgs doublet, with opposite hypercharge to that of the SM 
■ Higgs doublet, in order to cancel the gauge anomalies and generate masses for all the fermions of the Standard 
Model 3 . If we want broken supersymmetry to be effective in protecting the Higgs mass against large radiative 
corrections, then the supersymmetric partners of the Standard Model (SM) particles cannot be much heavier 



■ than about 1 TeV. The discovery of the superpartners of the SM particles is one of the main goals of present 
and future accelerators. 

When the SM gauge symmetry SU{2) x U{\) is broken, the fermionic partners of the two Higgs dou- 
blets {Hi, H2) of the MSSM mix with the fermionic partners of the gauge bosons, resulting in four neutralino 

■ states Xi J * = U 2, 3, 4, and two chargino states , j = 1, 2. In several models of low energy supersymmetry, 
H I the lightest neutralino is typically the lightest supersymmetric particle (LSP). In the MSSM, if we assume 

i?-parity (Rp) conservation, then the lightest supersymmetric particle is absolutely stable. There has been ex- 
tensive study of the neutralino states of the minimal supersymmetric standard model, and its extensions [3-[T2| , 
because the lightest neutralino, being the LSP, is the end product of any process that involves supersymmetric 
particles in the final state. In this work we will assume that the LSP is the lightest neutralino, it is stable, and 
that it escapes the collider experiments undetected. The composition and mass of the neutralinos and charginos 
will play a key role in the search for supersymmetry at high energy accelerators. The composition and mass 
of neutralinos will also determine the time-scale of their decays. The implications of mass patterns of the 
neutralinos in models with different particle content, or with specific SUSY breaking patterns were considered 
in some detail in [1, [ll|, [13] • 

At present an indirect phenomenological evidence for supersymmetry is obtained from the unification of the 
gauge couplings of the Standard Model in supersymmetric grand unified theories (GUTS) Furthermore, 
one of the most important prediction of grand unification is that of baryon number violating interactions, 
leading to proton decay. However, in supersymmetric grand unified theories proton decay is much slower than 
in nonsupersymmetric grand unified theories. The reason for this is that the unification scale in supersymmetric 
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GUTS is of the order of Mquts ~ fewlO"'^^ GeV, which is about 20 — 30 times larger than the corresponding 
scale in non-SUSY GUTS. Thus, proton decay via gauge boson exchange is negligible and main decay arises 
from dimension- 5 op erators with higgsino exchange. This leads to a rate of proton decay which is close to the 
existing bounds 'l5|. Furthermore, the range of neutrino masses as indicated by current experiments, when 
interpreted in terms of see-saw mechanism, point towards a large scale consistent with the scale Mquts of 
supersymmetric grand unification [l6| . 

It is, thus, natural to study the phenomenology of neutralino and charginos in an underlying grand unified 
theory. Most of the studies involving neutralinos and charginos in the minimal supersymmetric standard model 
have been performed with universal gaugino masses at the grand unification scale [17]. The neutralino and 
chargino masses depend on the soft SU{2)l and U{1)y gaugino masses M2 and Mi, the higgs(ino) parameter 
/Lt, and tan /3 = V2 /vi , where V2 and vi are the vacuum expectation values of the neutral components of the two 
Higgs doublets H2 and Hi. Most of the models assume the gaugino mass universality at the GUT scale, i.e. 
Ml = M2 = M3, where M3 is the SU{3)c soft gaugino mass. However, there is no specific theoretical reason 
for the choice of universal masses at the grand unification scale. It is possible to have nonuniversal gaugino 
masses in grand unified theories wherein the standard model gauge group is embedded in a grand unified gauge 
group. In a given supersymmetric model gaugino masses are generated from higher dimensional interaction 
terms involving gauginos and auxiliary parts of chiral superfields [l^. As an example, in SU{5) grand unified 
theory (GUT), the auxiliary part of a chiral superfield in higher dimensional terms can be in the representation 
1, 24, 75, or 200, or, in general, some combination of these representations. When the auxiliary field of one of 
the SU{5) nonsinglet chiral superfields obtains a vacuum expectation value (VEV), then the resulting gaugino 
masses are nonuniversal at the grand unification scale. Similar conclusions hold for other supersymmetric grand 
unified models. Furthermore, nonuniversal supersymmetry breaking masses are a generic feature in some of 
the realistic supersymmetric models. For example, in anomaly mediated supersymmetry breaking models the 
gaugino masses are not unified [H, [2^ , and hence are not universal. 

From above discussion it is clear, that the phenomenology of supersymmetric models depends crucially on 
the composition of neutralinos and charginos. This in turn depends on the soft gaugino mass parameters M2 
and Ml, besides the parameters and tan/?. Since most of the models discussed in the literature assume 
gaugino mass universality at the GUT scale, it is important to investigate the changes in the phenomenology 
of broken supersymmetry which results from the changes in the composition of neutralinos and charginos that 
may arise because of the changes in the pattern of soft gaugino masses at the grand unification scale [21]. The 
consequences of nonuniversal gaugino masses at the grand unified scale and the resulting change in boundary 
conditions has been considered in several papers. This includes the study of constraints arising from different 
experimental measurements (2^42^ . and in the study of supersymmetric dark matter candidates [25l [26j. 

One of the major goals of high energy colliders is to discover the supersymmetric partners of the Standard 
Model particles. In particular, a high energy e+e~ linear collider with a center-of-mass energy of y/s = 500 GeV 
in the first stage, and with a high luminosity C ~ 500 fb~\ will be important in determining the parameters 
of the broken sup ersy mmetric model with a high precision [27'-^3l'] . 

Recently in [I3, |32| a detailed study of the radiative production of neutralinos in electron-positron collisions 
in low energy supersymmetric models with universal gaugino masses at the grand unified scale has been carried 
out. In this paper we shall carry out a study of the implications of the nonuniversal gaugino masses, as they 
arise in grand unified theories, for the production of lightest neutralinos in electron-positron collisions. Since in 
a large class of models of supersymmetry the lightest neutralino is expected to be the lightest supersymmetric 
particle, it is one of the first states to be produced at the colliders. At an electron-positron collider, such as the 
International Linear Collider (ILC), the lightest neutralino can be directly produced in pairs [5, 33]. In collider 
experiments it will escape detection such that the direct production of the lightest neutralino pair is invisible. 

One must, therefore, look for the signature of neutralinos in electron-positron colliders in the radiative 
production process 

e+ + e-^x; + X?+7- (I-l) 

The signature of this process is a single high energy photon with missing energy carried away by the neutralinos. 
In this paper we carry out a detailed study of the process p. II) in supersymmetric grand unified theories with 
nonuniversal boundary conditions at the grand unified scale. Note that this process is suppressed by the square 
of the electromagnetic coupling. However, it might be the first process where the lightest supersymmetric 
particles could be observed at e+e" colliders. The process (jl.ip has been studied in detail in the minimal 
supersymmetric model [33 - l43 |. in various approximations. Calculations have also been carried out for MSSM 
using general neutralino mixing I42|-l44ll . This process has also been studied in detail in the next-to-minimal 
supersymmetric model (NMSSM) [l7ll32| . On the other hand different LEP collaborations (45l - l49j have studied 
the signature of radiative neutralino production in detail, but have found no deviations from the S M p rediction. 
Thus, they have only been able to set bounds on the masses of supersymmetric particles [isl - liTL l49|. 
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We recall that in the SM the radiative neutrino process 

e+e" ^ + i> + 7, (1.2) 

is the leading process with the same signature as p.ip . The cross section for the process p.2p depends on the 
number Ni^ of light neutrino species '50!| . This process acts as a main background to the radiative neutralino 
production process . There is also a supersymmetric background to the process coming from radiative 
sneutrino production 

e+e" ^ £> + £>*+ 7. (1.3) 
We shall consider both these processes as well, since they form the main background to the radiative process 

The layout of the paper is as follows. In Sec. |TT1 we review different patterns of gaugino masses that arise in 
grand unified theories. We will consider grand unified theories based on SU{b) and 5*0(10) gauge groups, and 
discuss the origin of nonuniversal gaugino masses in these models. Here we also calculate the elements of the 
mixing matrix which are relevant for obtaining the couplings of the lightest neutralino to the electron, selectron, 
and Z boson which control the radiative neutralino production process (11.11) . Furthermore, we also describe in 
detail the typical set of input parameters that are used in our numerical evaluation of cross sections. The set 
of parameters that we use are obtained by imposing various experimental and theoretical constraints discussed 
in Appendix |5] on the parameter space of the minimal supersymmetric standard model with underlying grand 
unification. These constraints will also be used throughout to arrive at the allowed parameter space for different 
models in this paper. 

Furthermore, we expect low energy observables from fiavor physics and 5 — 2 to have some impact in con- 
straining various models studied in this paper. The study of the impact of these observables on our analysis 
is, however, beyond the scope of present paper. 

In Sec. mil we summarize the phase space for the signal process, and also the cuts on the photon angle 
and energy that are used to regularise the infrared and collinear divergences in the tree level cross section. 
In Sec. IIVI we evaluate the cross section for the signal process (jl.l[) in different grand unified theories with 
nonuniversal gaugino masses, using the set of parameters obtained in Sec. |TT] for different patterns of gaugino 
mass parameters at the grand unified scale. We have included higher order QED radiative corrections in our 
calculations. We also compare and contrast the results so obtained with the corresponding cross section in the 
MSSM with universal gaugino masses at the grand unified scale. The dependence of the cross section on the 
parameters of the neutralino sector, and on the selectron masses is also studied numerically. 

In Sec.|V]we discuss the backgrounds to the radiative neutralino production process p.l[) from the SM and 
supersymmetric processes. An excess of photons from radiative neutralino production over the backgrounds 
measured through statistical significance is also discussed here and calculated for different grand unified models. 
Our results are summarized and the conclusions presented in Sec. IVII 

II. GAUGINO MASS PATTERNS IN GRAND UNIFIED THEORIES 

In this section we shall discuss soft supersymmetry breaking gaugino mass patterns that arise in SU{b) and 
S'O(IO) supersymmetric grand unified models, and the implications of these mass patterns for the neutralino 
masses and couplings. In Appendix [21 we summarize our notations for the neutralino mass matrix and the 
couplings of the lightest neutralino that are relevant to our studvflTl. [5ll|. Furthermore, in Appendix iBl we 
summarize the current experimental constraints (52| - [55j on the parameters of the neutralino mass matrix that 
we use in our calculations. 



A. Universal Gaugino Masses in Grand Unified Theories 

In supersymmetric models, with gravity mediated supersymmetry breaking, usually denoted as mSUGRA, 
the soft supersymmetry breaking gaugino mass parameters Mi, M2, and M3 and the respective gauge couplings 
Qi satisfy the renormalization group equations (RGBs) {\M^\ = Mg, the gluino mass) 

(f,l,-3), (II.l) 

(11.2) 



167r^ 



, dM,, 



dt 
dg, 
dt 



hot 
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at the one-loop order, where j = 1,2,3 refer to the U{1)y, SU{2)l and the SU{3) gauge groups, respectively. 
Here, gi = ^g', (72 = g, with g' and g as U{1)y, and SU{2)l gauge couplings, respectively, and (73 is the SU{Z)c 
gauge coupling. With the universal boundary conditions on the gaugino masses {ai = gfj'^i^, i — 1,2,3), we 
have 

Ml = M2 = A/3 = mi/2, (II.3) 
a\ = a2 = as = ac, (II-4) 

at the GUT scale Mq- Then the RGEs (|II.ip and (|II.2p imply that the soft gaugino masses scale like gauge 
couplings: 

Mi{Mz) ^ A'hjMz) ^ MsjMz) ^ rn^^ 
ai{Mz) ^ a2{Mz) ~ a^i^Mz) ~ aa ' 

The relation (III. 51) implies that out of three gaugino mass parameters only one is independent, which we are 
free to choose as the gluino mass Mg. The remaining soft gaugino mass parameters can be determined through 

ha 

A'h{Mz) = Tl^Mg ~ 0.14 M^, (II.6) 

3q;3 cos^ Ow 

A-hiMz) = Mg ^ 0.28 Mg, (II.7) 

az sm Ow 

where we have used the values of various couplings at the to be 

a'^{Mz) = m.Q, sin^ 61^1/ = 0.23, 03(^2) = 0.12. (II.8) 

For the gaugino mass parameters this leads to the ratio 

M 1 : M2 : Af3 ~ 1 : 2 : 7.1. (II.9) 

The gaugino mass parameters described above are the running masses evaluated at the electroweak scale Mz- 
Using the ratio (jll.9p and the constraint (|B.3|) . we have the lower bound on the parameter Mi 

Ml > 50 GeV, (11.10) 

in mSUGRA model. We shall implement this constraint on the parameter Mi in our calculations. 

For the case of universal gaugino masses at the grand unified scale, we shall use the set of parameters shown 
in TablelD We shall call this set of parameters the MSSM electroweak symmetry breaking scenario (EWSB) [56l |. 
This scenario has the advantage that it allows us to study the dependence of the neutralino masses and the 
radiative neutralino production cross section on /i and M2, and on the selectron masses. 

The values of different parameters in Table |T] have been arrived as follows. We first choose the smallest value 
of M3 to be around 1400 GeV as dictated by the experimental constraints on the gluino mass. We then vary 
it from 1400 to 3000 GeV in steps of 100 GeV. For the smallest value of M3 = 1400 GeV, we determine the 
values of soft gaugino mass parameters Mi and M2. Increasing the value of M3 in steps of 100 GeV, we obtain 
the corresponding values of Mi and M2 for different vales of M^ . This is shown in Fig. [T] as the curve labelled 
MSSM EWSB. Then we scan the values of the parameter Af2, corresponding to different values of M3, with 
values of the parameter fi varying between 110 to 200 GeV such that the mass of the lightest neutralino lies 
between 100 to 200 GeV. This is shown as contour plot in Fig. H For Table [J for the MSSM EWSB, we have 
chosen the values of M2 and /i in Fig. [5] which correspond to the lightest neutralino mass of 108 GeV and 
chargino mass larger than 110 GeV. Other values can be obtained by choosing a higher mass for the lightest 
neutralino. 

After selecting the values of the parameters for our analysis, we vary the value of tan /? in the range so 
that the top and bottom quark Yukawa couplings remain perturbative upto the grand unified scale. For this 
variation in the value of tan/3, we find that the mass of the lightest neutralino varies by only ± (2 — 3) GeV 
compared to its value for tan /? — 10, with very little change in the cross section for the radiative neutralino 
production. Since our analysis is insensitive to the value of tan /3, we have chosen tan /3 = 10 for definiteness. 

We note from the Appendix!^ the couplings of the lightest neutralino to electrons, selectrons, and Z bosons, 
which are used for the calculation of the radiative neutralino production cross section, are determined by the 
corresponding elements of the neutralino mixing matrix iVy . For the MSSM EWSB scenario of Table U the 
composition of the lightest neutralino is given by 

Nij = (0.345,-0.175,0.703,-0.596). (11.11) 
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FIG. 1: Dependence of the mass parameters Mi 
and M2 on the gluino mass Mg. Each triangle in 
the plot refers to a gluino of a particular mass, 
starting from 1400 GeV, and successively 
increasing by 100 GeV. 



FIG. 2: The lightest neutralino mass in fj, — M2 
plane in case of the MSSM EWSB scenario. The 

value of fi giving m^o « 108 GeV, for the 
smallest value of M2 satisfying the gluino mass 
constraint, is chosen for our analyses. 



tan /3 = 10 


= 130 GeV 


Ml = 197 GeV 


M2 = 395 GeV 


M3 = 1402 GeV 


At= 2300 GeV 


Ab= 2300 GeV 


A^= 1000 GeV 


m^o = 108 GeV 


m ± = 125 GeV 


mg^ = 156.2 GeV 


mp^ = 136 GeV 


m^o = 140 GeV 


m ± = 422.4 GeV 


mg^ = 156.7 GeV 


rrih = 124.3 GeV 



TABLE I: Input parameters and resulting masses of various states in MSSM EWSB scenario. 



The lightest neutralino in this scenario is dominantly a higgsino. 

From Table IXVI in Appendix \^ we see that for the composition of the neutralino in (III. lip , the neutralino - 
Z'^ coupling is enhanced compared to the coupling of the lightest neutralino with right and left selectrons Cu^l- 

As a benchmark, we shall use the radiative neutralino cross section for the MSSM with universal boundary 
condition (|II.3p for the gaugino mass parameters at the GUT scale. For this, we shall work with the set of 
parameters as shown in Table ID 



B. Nonuniversal Gaugino Masses in Grand Unified Theories 



We now consider the neutralino masses and mixing in the minimal supersymmetric standard model with 
nonuniversal boundary conditions at the GUT scale which arise in SU{5) and 50(10) grand unified theories. 
The masses and the composition of neutralinos and charginos are determined by the soft supersymmetry 
breaking gaugino masses Mi and M2, respectively, the supersymmetric Higgs (ino) mixing parameter /i, and 
the ratio of the vacuum expectation values of the two neutral Higgs bosons and H2, {H2) / {Hi) = tan/3. 
As discussed in subsection III Al in the simplest supersymmetric model with universal gaugino masses. Mi {i = 
1, 2, 3) are taken to be equal at the grand unified scale. However, in supersymmetric theories with an underlying 
grand unified gauge group, the gaugino masses need not necessarily be equal at the GUT scale. In this Section 
we briefly review the non universality of gaugino masses as it arises in SU{5) and 50(10) supersymmetric 
grand unified theories, and its implications for the neutralino masses and couplings. 

In grand unified supersymmetric models the soft gaugino masses are generated by coupling the field strength 
superfield W"" to /($), where /($) is the gauge kinetic function which is an analytic function of the chiral 
superflelds $ in the theory [l^ . The Lagrangian for the coupling of gauge kinetic function to the gauge field 
strength can be written as 

Cg.k. = j dHfabrnW^W" + h.c., (11.12) 
where a and b refer to gauge group indices, and repeated indices are summed over. The gauge kinetic function 
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fab{^) can be written as 

/ab($) = /o($^)'5a6 + 5]/„(<i>^)^ + (11.13) 

n 

Here $^ and the $" denote the singlet and the non-singlet chiral superfields, respectively. Also, /o(*i>*) and 
/„($'') are functions of gauge singlet superfields , and Mp denotes some large scale, e.g. the Planck scale. 
When the auxiliary part Fg, of a chiral superfield $ in /($) gets a VEV (Fg,), the interaction (|II.12[) gives rise 
to soft gaugino masses: 

Cg.k. ^ ^^A'^A^ + /i.e., (11.14) 

Alp 

where A"''' are gaugino fields. Here, we have denoted by A^, A^ and A'^ as the U{1), SU{2) and SU{3) gaugino 
fields, respectively. Since the gauginos belong to the adjoint representation of the gauge group, $ and F^ can 
belong to any of the representations appearing in the symmetric product of the two adjoint representations of 
corresponding gauge group. 

Since the SM can be embedded into a larger gauge group the question of unified gauge group needs to be 
discussed in order to study the implications for nonuniversal soft gaugino masses in grand unified theories. We 
recall that there is chain of group embeddings of the SM gauge group into a larger group [57| 

SUi3)c X SU{2)l X C/(l)y c SU{5) c SO{10) c Eq C Ej c Eg. (H.15) 

However, we note that in four-dimensional grand unified theories the gauge groups Ei and E^ do not lead 
to a chiral structure of the weak interactions, and hence are ruled out as grand unified gauge groups on 
phenomenological grounds. This leaves out only the three groups, SU{5), SO{10), and Eq as possible grand 
unified gauge groups in four dimensions. Here we shall study the implications of nonuniversal gaugino masses 
in the case of SU{5) and 5*0(10) grand unified gauge groups. 



1. SU(5) 



In this Section we shall consider the case of embedding of the SM gauge group into the grand unified gauge 
group SU{5). In the symmetric product of the two adjoint (24 dimensional) representations of SU{5), we have 



(24(g)24)sy„ 



1 ® 24 © 75 © 200. 



(H.16) 



In the simplest case where <& and Fg, are assumed to be in the singlet representation of SU{5), we have equal 
gaugino masses at the GUT scale. But, as we can see from (jII.16|) . $ and Eg, can belong to any of the non-singlet 
representations 24, 75, and 200 of SU{5). In such cases the so ft g augino masses are unequal but related to 
one another via the representation invariants of the gauge group [2l|. In Table |TT] we show the ratios of gaugino 
masses which result when Eg, belongs to different representations of SU{5) in the decomposition (jII.16p . In 
this paper, for definiteness, we shall study the case of each representation independently, although an arbitrary 
combination of these is also allowed. 

In the one-loop approximation, the solution of renormalization group (RG) equations for the soft supersym- 
metry breaking gaugino masses Mi, M2, and M3 can be written as [58^ 



M,(t) ^ M,(GUT) 
a^{t) a,(GUT) ' 



1,2,3. 



(11.17) 



Then at any arbitrary scale we have : 



3 cos^ 



Ml (GUT) \ 
ai(GUT)y' 



M, = 



a /M2(GUT)\ 



sin'^ew Va2(GUT) ; 



M3 = as 



M3(GUT) \ 
as (GUT) J 



(H.18) 



With these results, we can write the gaugino masses for the 24 dimensional representation of SU{5) as 



7 



SU{5) 


Mi- 












1 


1 


1 


1 


1 


2 


7.1 


24 


1 


3 


-2 


1 


6 


-14.3 


75 


1 


3 
5 


1 

5 


1 


-1.18 


-1.41 


200 


1 


1 


1 

in 


1 


0.4 


0.71 



TABLE II: Ratios of the gaugino masses at the GUT scale in the normaUzation Mi{GUT) — 1, and at the 
electroweak scale in the normalization Mi{EW) — 1 for F-terms in different representations of SU{5). These 
results are obtained by using 1-loop renormalization group equations. 



Similarly, for the 75 dimensional representation of SU{5), we have the result 

Ml ^ a \ / 1 \ M2 ^ / a \fl 

Ms \3 cos^ Ow J \a3 J ' M3 \smHw J K^s 

and finally for the 200 dimensional representation of SU{5) we have 

- - 10 (-^] f-V ^ = 2 (^] (^) . (11.21) 

To compute the ratios of the gaugino masses at the electroweak (EW) scale Mz for different representations of 
SU{b) in the product pLT6| . we use the relevant renormalization group equations for the soft gaugino masses. 
In Table In] we show the approximate results for these masses at the electroweak scale Mi{EW). These are 
calculated using one loop renormalization group equations for the gaugino masses and the gauge couplings. 
The effect of two- loop calculations is to increase the ratio M1/M2 by an amount which is not significant. It is 
important to note that these results are consistent with the unification of gauge couplings 

af = af ^af a^i^ 1/25), (11.22) 

at the GUT scale, where we have neglected the contribution of non universality in gaugino masses to the gauge 
couplings, which is not significant. 



(11.20) 




FIG. 3: The lightest neutralino mass in fi — M2 
plane for SU (5) with $ and F$ in the 24 
dimensional representation. The value of ^ 
giving m^o « 108 GeV, for the smallest value of 
M2 satisfying the gluino mass constraint, is 
chosen for our analyses. 



FIG. 4: The lightest neutralino mass in — M2 
plane for SU{5) with $ and ^$ in the 75 
dimensional representation. The value of fi 
giving m^o « 108 GeV, for the smallest value of 
M2 satisfying the gluino mass constraint, is 
chosen for our analyses. 



The input parameters for 24, 75 and 200 dimensional representations of SU{5) in Table HIl are chosen in a 
manner similar to the one we used for the parameter space of MSSM EWSB. The lowest values of Mi and M2 
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FIG. 5: The lightest neutrahno mass in fi ~ M2 plane for SU{5) with $ and i^* in the 200 dimensional 
representation. The value of /i giving m^o ss 108 GeV, for the smallest value of AI2 satisfying the gluino mass 

constraint, is chosen for our analyses. 



tan /3 = 10 


= 138 GeV 


Ml = 149 GeV 


Ma = 890 GeV 


M-i = -2121 GeV 


At= -1000 GeV 


Ab= -2050 GeV 


A^^ -2050 GeV 


m^o = 108 GeV 


m^± = 138.5 GeV 


mg^ = 156 GeV 


mp^ = 136 GeV 


m^o = 146 GeV 


m ± = 903 GeV 


= 157 GeV 


mh = 123.5 GeV 



TABLE III: Input parameters and resulting masses for various states in SU{5) supersymmetric grand unified 
theory with $ and _Flj, in the 24 dimensional representation. 
We shall refer to this model as [5'f/(5)]24 in the text. 



satisfying the gluino mass constraint are chosen from Fig. [T] corresponding to various SU{5) representations 
shown there. We then plot the values of M2 versus fi for various representations in contour plots shown in 
Figs. [3l H] and [5l The values of M2 and fj, are then selected which result in the lightest neutralino mass of 108 
GeV and chargino mass larger than 110 GeV. For SU{5)r5 and SU{5)2qq the bino and wino mass parameters 
Ml, M2 ^ fi, therefore the LSP mass in these cases is almost equal to fi. This is seen from Figs. H] and [5l where 
the contours in the M2 ~ fi plane are independent of M2 ■ 

The input parameters and the resulting masses for the 24, 75, and 200 dimensional representations of SU{5) 
which result in nonuniversal gaugino masses at the grand unified scale obtained in a manner described above 
are shown in Tables Hill IIVI and fVl respectively. In arriving at the parameter values in these Tables, we have 
taken into account various theoretical and phenomenological constraints, including the electroweak symmetry 
breaking at the correct scale, as described in the Appendix |B] Other values can be obtained by choosing larger 
values of the parameter M3. 

The composition of the lightest neutralino for the different representations of SU{5) in Table HIl is obtained 
from the mixing matrix for the choices of parameters given in Tables lllli IIVI and [V] This composition is 
calculated to be: 

1. SU{5) with $ and F$ in the 24 dimensional representation (labelled as model [SU{5)]24)- 

Nij ^ (0.638, - 0.054, 0.598, - 0.482); (11.23) 



tan /3 = 10 


^l = 108 GeV 


Ml = 993 GeV 


Ma = -1172 GeV 


M3 = -1400 GeV 


At^ 1000 GeV 


Ab^ 2300 GeV 


A.r= 2300 GeV 


m^o = 108 GeV 


m ± = 109 GeV 


mg^ = 156 GeV 


mc^ = 136 GeV 


m^o = 110 GeV 

X2 


m ± = 1916 GeV 

X, 


ma^ = 157 GeV 


ruh = 123.5 GeV 



TABLE IV: Input parameters and resulting masses for various states in SU{5) supersymmetric grand unified 
theory with $ and F$ in the 75 dimensional representation. 
We shall refer to this model as [SU{5)]75 in the text. 
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tan/3 = 10 


/i = 111 GeV 


Ml = 1970 GeV 


Ma = 788 GeV 


M-i = 1399 GeV 


At= 1000 GeV 


2300 GeV 


yl,= 2300 GeV 


m^o = 107.7 GeV 


m ± = 111 GeV 


mg^ = 166 GeV 


mc^ = 136 GeV 


m^o = 117 GeV 


m ± = 807 GeV 

X, 


mg^ = 157 GeV 


ruh = 123.7 GeV 



TABLE V: Input parameters and resulting masses for various states in SU (5) supersymmetric grand unified 
theory with $ and F$ in the 200 dimensional representation. 
We shall refer to this model as [S'I7(5)]200 in the text. 

2. S'[/(5) with <i> and _F$ in the 75 dimensional representation (labelled as model [SU {h)]^^^): 

Nij = (0.031, 0.056, - 0.710, -0.700); (11.24) 

3. SU{b) with and in the 200 dimensional representation (labelled as model [S'J7(5)]2oo): 

Nij = (0.018, - 0.085, 0.719, - 0.689). (11.25) 

We note from pi.23|) . (jll.24[) . and (IIL25P that for the 24 dimensional representation of S'(7(5), the dominant 
component of the ncutralino is the bino, whereas for the other representations of SU{5), there is a higgsino like 
lightest neutralino. Thus, for 75 and 200 dimensional representations the neutralino being higgsino like couples 
weakly to the selectron, with the dominant contribution to the cross section coming from the neutralino- Z° 
coupling. 

2. 50(10) 

We now consider the embedding of the SM gauge group in a 5*0(10) supersymmetric grand unified theory. 
Since the adjoint representation of 50(10) is 45 dimensional, $ and Fg, can belong to any of the following 
representations appearing [59^ in the symmetric product of two 45 dimensional representations of 50(10): 

(45 X 45)symm = 1 ® 54 © 210 ® 770. (11.26) 

There are three maximal proper subgroups of 5*0(10) which are consistent with the fermion content of the 
Standard Model. These are (i) 5C/(5) C 50(10) with the normal (nonfiipped) embedding; (ii) SU{5y x U{1) C 
5*0(10 ) wit h the f lipped embedding; and {in) 5*1/(4) x 5*C/(2)l x SU{2)r C 5*0(10) embedding. Using relations 
(|IL17p and (|II.18p . we obtain the gaugino mass parameters at the GUT scale for different representations that 
arise in the symmetric product of two adjoint representations of 50(10) with the relevant embedding of the 
SM gauge group in 50(10). These are shown in Tables IVI[ IVIIi and I Villi The ratio of the gaugino masses at 
the GUT scale for the different cases for 50(10) shown in these Tables can be scaled down to the electroweak 
scale, as described above. The results for the gaugino masses at the electroweak scale for these cases are also 
shown in Tables ED EIIl andEUH 

We note from Table IVll that the ratios of gaugino masses for the different representations of 5*0(10) in the 
symmetric product pi.26|) with the unflipped embedding 5'C/(5) C 5*0(10) are identical to the corresponding 
gaugino mass ratios in Table [TT] for the embedding of SM in SU{h). Therefore, the input parameters and 
the resulting masses for the gaugino mass ratios in Table I VII for 50(10) are identical to the corresponding 
Tables ImlllVl andElfor 5*f/(5). 

On the other hand, for the flipped embedding 5*C/(5)' x [/(I) C 50(10), Table IVlIl the gaugino mass ratios 
for the 210 and 770 dimensional representations of the grand unified gauge group can be different from the 
corresponding ratios for 5*[/(5). 

In the case of 5*0(10 we choose the parameter values in a manner similar to the case of MSSM EWSB and 
the 5*t/(5) grand unified theory. The result of such a procedure is shown in Fig. [l] and in Figs HI [7] and HI The 
resulting values of parameters are shown in Tables IIXI 1X1 and IXII 

We note that the ratio of the gaugino masses for the (24, 0), (75, 0) of the 210 dimensional representation of 
5*0(10), and (75, 0), (200, 0) of the 770 representation of 50(10) follow the same pattern as the (24, 0) of the 
54 dimensional representation of 50(10). The behaviour of the cross section for the radiative neutralino cross 
section will be similar in these cases due to the fact that the dominant component of the lightest neutralino in 
this case is a bino. Similarly, the (24, 0) of the 770 dimensional representation of 5*0(10) has the same pattern 
as 200 dimensional representation of SU{5). Here also the behavior will be same. Because of this we shall 
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^□(lO) 


SU{5) 








Mf 




Mi"""' 


1 


1 


1 


1 


1 


1 


2 


7.1 


54 


24 


1 


3 


-2 


1 


6 


-14.3 


210 


1 


1 


1 


1 


1 


2 


7.1 




24 


1 


3 


-2 


1 


6 


-14.3 




75 


1 


3 
5 


1 

5 


1 


-1.18 


-1.41 


770 


1 


1 


1 


1 


1 


2 


7.1 




24 


1 


3 


-2 


1 


6 


-14.3 




75 


1 


_ 3 
5 


1 

5 


1 


-1.18 


-1.14 




200 


1 


1 
5 


1 

in 


1 


0.4 


0.71 



TABLE VI: Ratios of the gaugino masses at the 
GUT scale in the normahzation Mi{GUT) = 1, 
and at the electroweak scale in the normalization 
Ml (EW) = 1 for F-terms in representations of 
SU{5) C SO{10) with the normal (nonflipped) 
embedding. These results have been obtained at 
the 1-loop level. 



SO(IO) 


[SU{5y X U{l)]fHpped 


Mt 












1 


(1,0) 


1 


1 


1 


1 


2 


7.1 


54 


(24,0) 


1 


3 


-2 


1 


6 


-14.3 


210 


(1,0) 
(24,0) 


1 
1 


5 

19 
15 
7 


5 

19 
10 
7 


1 
1 


-0.52 
-4.2 


-1.85 
10 




(75,0) 


1 


-15 


-5 


1 


-28 


-33.33 


770 


(1,0) 
(24,0) 
(75,0) 


1 
1 
1 


5 

77 
15 
101 

-15 


5 

77 
10 
101 

-5 


1 
1 
1 


0.13 
0.3 

-28 


0.46 
-0.70 
-33.3 




(200,0) 


1 


5 


5 
2 


1 


9.33 


16.67 



TABLE VIL Ratios of the gaugino masses at the 
GUT scale in the normalization Mi{GUT) = 1, and 
at the electroweak scale in the normalization 
Mi{EW) = 1 at the 1-loop level for F-terms in 
representations of flipped SU{^)' x U{1) C SO{lQ). 



50(10) 


SU{A) X SU{2)r 












M3^^ 


1 


(1,1) 


1 


1 


1 


1 


2 


7.1 


54 


(1,1) 


1 


3 


2 


1 


6 


-14.3 


210 


(1,1) 


1 


5 
3 





1 


-3.35 







(15,1) 


1 





5 
4 


1 





-9.09 




(15,3) 


1 








1 








770 


(1,1) 


1 


25 
19 


10 
19 


1 


2.6 


3.7 




(1,5) 


1 








1 










(15,3) 


1 








1 










(84,1) 


1 





5 


1 





1.11 



TABLE VIII: Ratios of the gaugino masses at the GUT 
scale in the normalization Mi{GUT) = 1, and at the 
electroweak scale in the normalization Mi (EW) = 1 at the 
1-loop level for i^-terms in representations of 
SU{4) X SU{2)l X SU{2)r c 50(10). 



tan/3 = 10 


M = 116 GeV 


Ml = 760 GeV 


Ma = -395 GeV 


Ms = -1405 GeV 


At= 1000 GeV 


At= 2300 GeV 


Ar^ 2300 GeV 


m^o = 108 GeV 


m ± = 111 GeV 


me„ = 156 GeV 


mp^ = 136 GeV 


m^o = 122 GeV 


m ± = 422 GeV 


mg^ = 158 GeV 


rrih = 124.2 GeV 



TABLE IX: Input parameters and resulting masses for various states in SU{5Y x U{1) C 50(10) 
supersymmetric grand unified theory with $ and in the 210 dimensional representation with 
5C/(5)' X U{1) in (1,0) dimensional representation. 
We shall refer to this model as [5'0(10)]2io in the text. 



focus on the (1, 0) option for the 210 and 770 dimensional representations in Table rvlll Similarly, because of 
phenomenological reasons, we focus on the (1, 1) of SU{A) x SU{2)fj of the 770 representation only. 

For the parameters of Tables HXl 1x1 and IXII the composition of the lightest neutralino is given by the following. 

1. 50(10) where 5[/(5)' x [/(I) c 50(10) with $ and in the 210 dimensional representation with 
5t/(5)' X U{1) in (1,0) dimensional representation (labelled as model [50(10)]2io): 

Nij = (0.038, 0.194, -0.719, - 0.666); (11.27) 



2. 50(10) where 5C/(5)' x U{1) C 50(10) with $ and F^ in the 770 dimensional representation with 
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800 



700 




120 140 160 180 200 



fi (GeV) 

FIG. 6: The ligh test nc utralino mass in ^ - Ma plane for S'O(IO) where SU{b)' x U{1) C S'O(IO) with $ and 
in the 210 ( III. 271) . The value of /i giving m^o ~ 108 GeV, for the selected value of M2 is considered for 

our analyses. 



800 
700 



> 




120 140 160 180 200 

// (GeV) 



FIG. 7: The lightest neutralino mass in — M2 
plane for 6*0(10)770 scenario. The value of ji 

giving m^o w 108 GeV, for the smallest value of 
M2 satisfying the gluino mass constraint, is 
chosen for our analyses. 



2000 - 

1800 
> 1600 




120 140 160 180 200 

1^ (GeV) 



FIG. 8: The lightest neutralino mass in — M2 
plane for 5*0(10)770' scenario. The value of /i 
giving TO^o sa 108 GeV, for the smallest value of 
M2 satisfying the gluino mass constraint, is 
chosen for our analyses. 



SU{5y X U{1) in (1,0) dimensional representation (labelled as model [S'0(10)]77o): 

Nij = (0.011, - 0.193, 0.724, - 0.663); (11.28) 

3. 5*0(10) where 5C/(4) x SU{2)r x SU{2)l C 50(10) with $ and F$ in the 770 dimensional representation 
with 5?7(4) X SU{2)r in (1,1) dimensional representation (labelled as model [5'O(10)]770'): 

Nij = (0.125, - 0.660, 0.721, - 0.678), (11.29) 



tan /3 = 10 


^ = 118 GeV 


Ml = 3038 GeV 


Ma = 395 GeV 


M3 = 1398 GeV 


At= 1000 GeV 


Ab= 2300 GeV 


Ar^ 2300 GeV 


m^o = 108 GeV 


m ± = 113 GeV 


mg„ = 156 GeV 


mc^ = 136 GeV 


m^o = 125 GeV 


m ± = 422 GeV 


mg^ = 157 GeV 


mh = 124 GeV 



TABLE X: Input parameters and resulting masses for various states in SU{5y x U{1) C 5*0(10) 
supersymmetric grand unified theory with $ and in the 770 dimensional representation with 
5[/(5)' X U(l) in (1,0) dimensional representation. 
We shall refer to this model as [5'0(10)]77o in the text. 
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tan /3 = 10 


/i = 113 GeV 


Ml = 378 GeV 


Mi = 985 GeV 


M-i = 1402 GeV 


At^ 1000 GoV 


2300 GeV 


2300 GeV 


m^o = 107.4 GeV 


m ± = 114 GeV 


TTieH = 156 GeV 


mf, = 136 GeV 


m^o = 120 GeV 


m ± = 1000 GeV 


mg^ = 157 GeV 


nih = 123.5 GeV 



TABLE XI: Input parameters and resulting masses for varfous states in SU{4) x SU{2)]i x SU{2)l C 5*0(10) 
supersymmetric grand unified theory with $ and in the 770 dimensional representation with 
SU{A) X SU{2)ji in (1,1) dimensional representation. 
We shall refer to this model as [S'O(10)]770' in the text. 

implying thereby that higgsino is the dominant component for the 210 and 770 dimensional representations 
with the embedding SU{b)' x U{1) C 5*0(10) and for the 770 dimensional representation with the embedding 
5*1/(4) X SU(2)r X SU{2)l C 50(10). From the Figs El [7] and [8] we see that, because the lightest neutralino is 
dominantly a higgsino, the contours are almost independent of M2. 

Thus, in these cases the dominant contribution to the radiative neutralino production cross section will come 
from the neutralino-Z° coupling. Since the LSP for most of the scenarios considered here has a dominant 
higgsino component, the Z° width imposes a strict constraint, as the Z° decay rate involves coupling to the 
higgsino component of the neutralino. We have imposed the LEP constraint on the anomalous Z° decay width 
in our calculations : 

r(Z ^ xlxl) < 3MeV. (11.30) 

III. RADIATIVE NEUTRALINO PRODUCTION IN GRAND UNIFIED THEORIES 

In this Section We shall calculate the cross section for the radiative neutralino production process 

e-(pi) + e+{p2) ^ X?(fci) + X?(fc2) + 7(9), (ni.l) 

in SU{5) and 50(10) grand unified theories with nonuniversal gaugino masses at the grand unified scale. The 
symbols in the brackets denote the four momenta of the corresponding particles. At the tree level, the Feynman 
diagrams contributing to the radiative neutralino production process are shown in Fig. ^ In order to calculate 
the cross section for the radiative process pil.ll) . we require the couplings of the neutralinos to electrons, the 
selectrons, and to the bosons. These couplings are obtained from the neutralino mixing matrix (IA.2P as 
in the Appendix \^ with the values of the soft SUSY breaking gaugino mass parameters Mi and M2 for the 
respective grand unified theory, as calculated in Section HIl 

The elements of the neutralino mixing matrix Nij for the case of 5*C/(5) and 5*0(10), which are relevant 
for our calculations, were calculated in the previous section. As indicated earlier, as a benchmark, we shall 
calculate the radiative neutralino cross section for the MSSM with universal boundary condition pi. 31) for 
the gaugino mass parameters at the GUT scale, for which we will work with the parameters in the MSSM 
electroweak symmetry breaking scenario (EWSB) [s^. This set of parameters is summarized in Table HI The 
relevant elements of the neutralino mixing matrix Nij are summarized in pi. Ill) . 

A. Cross Section for the Signal Process 

From Fig. [9l we see that the t- and w-channel exchange of right and left selectrons eu.L, and via Z boson 
exchange in the s channel contribute to the process pil.ip . The differential cross section for the process pil.ip 
can be written as [H, HO] 

where p/ and Ef are the final three-momenta ki, k2, q and the final energies i^xu ^X2^ and of the 
neutralinos and the photon, respectively. Using the standard technique, we average over initial spins and sum 
over the spins of the outgoing neutralinos. We also sum over the polarizations of the outgoing photon. Putting 
all this into effect, the squared matrix element |A^p in pil.2[) can be written as [1^ 



(III.3) 
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FIG. 9: Feynman diagrams contributing to the radiative neutralino production e+e~ — XiXi7- There are six 
other diagrams which are exchange diagrams corresponding to (i, ii, in, vi, vii, viii), with w-channel 
exchange of selectrons, wherein the neutrahnos are crossed in the final state. 

where Tij are squared amphtudes corresponding to the Feynman diagrams in Fig. [S] The phase space for the 
radiative neutrahno production process in (IIII.2I) is described in detail in (36l |. 



B. Radiative Corrections 



The next generation linear colliders are designed to have a high luminosity, which in turn will require beams 
with bunches of high densities. Due to the above requirement of high number density in the bunch, there 
arise problems with the generation of strong electromagnetic fields in and around the bunches. This in turn 
generates initial state radiation (ISR) and beamstrahlung effects. These effects have been studied extensively 
in the past [6ll - l64| in the context of the future linear colliders. Therefore ISR and beamstrahlung effects have 
to be considered in any realistic calculation of the cross sections at a future linear collider since they result in 
the loss of beam energy along with the disturbance of the initial beam calibration. 

The majority of the emitted photons are soft and are lost down the beam pipe. Only the hard photons can 
be tagged. The radiated hard photon usually carries from the radiating particle a certain amount of energy, 
resulting in the energy distribution of the initial beams. So a precise knowledge of ISR along with distribution 
of the photon spectrum resulting from beamstrahlung and the behaviour of the electrons, positrons after the 
emission is required. We have calculated the radiative effects for our process and the background processes 
using CalcHEP 56] , with parameters given in Table IXIII |65l] . The resulting spectrum of electrons has been 
calculated using the structure function of the initial leptons valid upto all orders in perturbation theory. We 
note that the radiative effects are included in all our calculations of the signal and background processes. 



IV. NUMERICAL RESULTS 



The tree- level cross section for radiative neutralino production (jlll.ip . and the standard model background 
from radiative neutrino and sneutrino production, (jL2]) and (jOj, have been calculated using the program 
CalcHEP [5^ . As noted above we have included the effects of radiative corrections to the signal as well as the 



14 



Collider Parameters 


ILC 


<^x (nm) 


640 


ay (nm) 


5.7 




300 


N (10^0) 


2 



TABLE XII: Beam parameters for the ILC, where N is the mimber of particles in the bunch and Ux, Uy are 
the transverse bunch sizes at the interaction point, with Uz as the bunch length. 





aEWSB 


aEWSB + R.E. 


O" SM 


asM + R-E. 


GeV 


(fb)xlO-" 


(fb)xlO"' 


(fb)xlO^ 


(fb)xlO^ 


300 


0.9072 


0.7586 


2.1899 


2.4187 


400 


1.4035 


1.2821 


2.3691 


2.4373 


500 


1.3963 


1.3377 


2.4191 


2.4321 


600 


1.2540 


1.2383 


2.4329 


2.4518 


700 


1.0922 


1.1030 


2.4226 


2.4359 


800 


0.9443 


0.9721 


2.3980 


2.3648 


900 


0.8179 


0.8568 


2.3687 


2.3318 


1000 


0.7124 


0.7576 


2.3341 


2.2934 



TABLE XIII: Cross section of the signal process in MSSM EWSB scenario, along with the the SM irreducible 
background with and without the inclusion of the radiative effects (R.E.) due to ISR and beamstrahlung. 

background processes. We recall that the tree level cross sections have infrared and collinear divergences, which 
must be regularized by imposing cuts on the fraction of beam energy carried by the photon and the scattering 
angle of the photon 36] . To implement this regularization, we define the fraction of the beam energy carried 
by the photon as a; = i?7/-Ebcam, where = 2ii'bcam is the center of mass energy, and is the energy carried 
away by the photon. We then impose the following cuts on x, and on the scattering angle 9-^ of the photon [66| : 

rn?a 

0.02 < a; < 1 - -^r^. (IV.l) 

beam 

-0.95 < cosO-y < 0.95. (IV.2) 

Note that the lower cut on x in pV.ip corresponds to photon energy Ej = 5 GeV for = 500 GeV. The 
upper bound of (1 — 771^0 /-^beam) *^ corrcspouds to the kinematical limit of radiative neutralino production 
process. The detector acceptance cut on the photon is applied so as to enhance the signal over the main 
irreducible background coming from the radiative neutrino production. We did not find any other cuts which 
would significantly reduce the background. 

The mass of the lightest neutralino for MSSM, with universal boundary conditions on the gaugino mass 
parameters at the GUT scale is taken to be m^o = 108 GeV from the EWSB scenario, to implement the cuts 
on the photon energy in the calculation of the cross sections. Using Eq. (jIV.l|) we get a fixed upper limit 
^max ^ 203.3 GeV for MSSM at = 500 GeV for the photon energy. We have used these cuts for both signal 
and background processes for all the scenarios that we have considered in this paper. 
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FIG. 10: The photon energy distribution da/dEj 
for the radiative neutrahno production including 
radiative effects for the SU{5) grand unified 

theory with nonuniversal gaugino masses 
compared with MSSM EWSB with universal 
gaugino masses. 



FIG. 11: The photon energy distribution da/dE. 

including radiative effects for the radiative 
neutralino production for 5*0(10) grand unified 
theory with nonuniversal gaugino masses 
compared with MSSM EWSB with universal 
gaugino masses. 



A. Photon Energy (E^) Distribution and Total Beam Energy (y^) Dependence 



To begin with, we have calculated the energy distribution of the photons from radiative neutralino production 
in MSSM EWSB scenario as well as for the different GUT models considered in this paper. 
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FIG. 12: Total cross section a for the radiative 
process e+e" XiXi7) with the inclusion of 
radiative effects as a function of center of mass 
energy for the SU (5) grand unified theory 
with nonuniversal gaugino masses compared with 
MSSM EWSB scenario with universal gaugino 
masses at the grand unified scale. 



FIG. 13: Total cross section a for the radiative 
process e+e~ — >■ XiXi7j with the inclusion of 
radiative effects as a function of center of mass 
energy ^/s for the 50(10) grand unified theory 
with nonuniversal gaugino masses compared with 
MSSM EWSB scenario with universal gaugino 
masses at the grand unified scale. 



In Table IXlIIl we illustrate the changes in the signal cross section with and without the inclusion of radiative 
effects for the MSSM EWSB scenario, along with the irreducible SM background with the application of cuts 
described earlier. The radiative effect leading to the distribution of the beam energy can be distinctly seen from 
this Table, with the increase of cross section at higher cm. energies and the corresponding decrease at lower 
cm. energies for the signal process. Similar behaviour holds in other models for the signal process and the 
supersymmetric background processes. In case of the background process, the radiative neutrino production, 
the presence of the Z pole along with the massless particles in the final state leads to a decrease of cross section 
at higher cm. energies and vice versa. The plots here only show the results with the inclusion of radiative 
effects of ISR and beamstrahlung, since the behaviour in case of different scenarios for the signal process and 
supersymmetric background with and without the inclusion of radiative effects is similar to the MSSM EWSB 
scenario shown in Table IXlIIl 

In Figs. [To] and [TT] we show the energy distribution of photons for models with nonuniversal gaugino masses 
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in grand unified theories based on SU{5) and SO (10), where we have included the radiative effects. We have 
compared this with the photon energy distribution for the MSSM EWSB model with universal gaugino masses 
at the GUT scale. The energy dependence of the total cross section for these models is also calculated and 
this is shown in Figs. [T2] and [T3l From Figs. [TOl and [Til as well as from Figs. [T2] and [T3l it can be seen that the 
signal in case of MSSM EWSB and [SU{5)]24 is enhanced compared to the other scenarios considered here. 
This is mainly due to the fact that the dominant component of the neutralino in [SU{5)]24 is a bino, whereas 
in other cases the lightest neutralino is dominantly a higgsino state. The MSSM EWSB scenario predicts a 
lightest neutralino with a dominant higgsino component, but it also has a significant bino component leading 
to the enhancement of right selectron-electron-neutralino coupling. For the other cases with a higgsino like 
neutralino the t- and u- channel exchange of e^.L, is suppressed, with the only contribution coming from off 
shell Z decay. In order to find ways of enhancing the signal, one must study the dependence of the signal on 
selectron ejix masses as well as on the parameters fi and M2, which determine the neutralino mixing elements. 



B. Dependence on ^ and M2 

As can be seen from the neutralino mass matrix in the Appendix |BJ the mass of the lightest neutralino 
depends on the parameters fi and M2. Therefore, it is important to study the dependence of cross section 
of the signal process on these parameters. Since fi and M2 are independent parameters, we have studied the 
dependence of the cross section a{e^e~ — ^ XiXi7) on these parameters independently. We have considered 
here all the scenarios with both universal and nonuniversal gaugino masses. The values of the parameters fi and 
M2 are chosen so as to avoid color and charge breaking minima, unbounded from below constraint on scalar 
potential, and to satisfy phenomenological constraints on different sparticle masses discussed in Appendix iBl 

We note here that we have carried out a check on the parameter space used in our calculations as to whether 
the complete scalar potential has charge and color breaking minima (CCB) which is lower than the electroweak 
minimum. We have also checked whether the scalar potential is unbounded from below (UFB). The criteria 
used for these conditions are 

Aj < 3(m^^ + {mj^ +ti'^+mH,), (IV. 3) 

mH,+mH, > 2\B^^\, (IV.4) 

respectively, at a scale > M|-yygg. Here / denotes the fermion generation, and A is the trilinear supersym- 
metry breaking parameter. We have implemented these conditions through the SuSpect package [67[ which 
computes the masses and couplings of the supersymmetric partners of the SM particles. For each model consid- 
ered in this paper, we perform the renormalization group evolution to calculate the particle spectrum. While 
doing so we check for the consistency of the chosen parameter set with electroweak symmetry breaking and 
that the conditions (lIV.Sp and (jIV.4p are satisfied. 
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FIG. 14: The total radiative neutralino 
production cross section a with radiative effects 
included as a function of /i in the range /i € [110, 
160] GeV for different models considered in this 
paper at ,/s = 500 GeV. 
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FIG. 17: Total cross section a along with 
radiative effects for the radiative neutralino 
production versus at y/s = 500 GeV. 



In Fig. [T3] we show the /i dependence of the cross section for different models considered in this paper. For 
a wide range of values of /x, the [5'C/(5)]24 scenario satisfies all experimental constraints, with x° as the LSP. 
Since the neutralino in this case is mainly a bino like state, with values of /i > Mi, the neutralino mass is 
relatively insensitive to the values of /i. For the other scenarios with a higgsino type lightest neutralino, the 
cross section is sensitive to the value of /i. Since m-j^o oc /i, above a certain value of /i, Xi ceases to be the 
lightest supersymmetric particle. Depending on the percentage of the higgsino component, the cross section 
changes with value of ji. In [S'J7(5)]75, since Afi,Af2 ^ w^? depends on ^, so from Fig. [T31 it can be seen 
that the signal cross section decreases with /i due to increasing m^o . Most of the scenarios considered here 
are tightly constrained as a function of /x, with neutralino as the LSP. This is due to the various limits on the 
sparticles masses from the experiments. The cross section for some scenarios in this region is too small to be 
observed at ILC with ^/s — 500 GeV, even with an integrated luminosity of 500 fb~^. 

We next show in Fig. [15] the dependence of the radiative neutralino cross section on the soft gaugino mass 
parameter M2 for different models that we have studied in this paper. From this Fig. we note that the total 
cross section decreases with increasing values of M2. The elements of the neutralino mixing matrix changes 
with the variation of the wino parameter M2, which in turn change the values of the vertices contributing to 
the radiative neutralino production, Table IXVl in the Appendix El The range of M2 considered for the different 
models is taken from Fig. [1] The variation of the cross section with M2 for S'[/(5)75 and 50(10)770' are shown 
in the inset of Fig. [151 their M2 range is different from the other scenarios. A lower value of M2 favours a 
cross section which can be measured experimentally for some scenarios, whereas for higher values the difficulty 
in measurement of cross section increases. 



C. Dependence on selectron masses 

The signal process cr(e"'"e~ XiXi7) mainly proceeds via right and left selectron cr^l exchange in the t- and 
M-channels. Different models that we have considered in this work have the selectron masses as independent 
parameters. Figs. [1^1 and [T71 show the dependence of the total cross section a for the radiative neutralino 
production on the left and right selectron masses. From Fig. [THl we see that the cross section is not sensitive 
to the left selectron mass for [SU{b)]2i scenario, because the LSP is predominantly a bino. As a consequence, 
the dominant contribution to the cross section comes from the right selectron exchange, with the contribution 
from the left selectron exchange suppressed by tan 0w in the coupling, as can be seen from Table IXVl in the 
Appendix [21 As a consequence the cross section has very little sensitivity to the left selectron mass, but is 
sensitive to the right selectron mass in the range 150 - 1000 GeV. The MSSM EWSB scenario has a LSP with 
a dominant contribution coming from higgsino, but has a significant bino contribution as well. Due to this the 
cross section is relatively insensitive to mass of cl, but a sensitivity is seen with respect to right selectron mass. 

In the other scenarios with a higgsino type LSP, the neutralino coupling with cr.l is suppressed, making 
these almost insensitive to the selectron masses. As can be seen from the Figs. [TBI and ITTI the cross section for 
these scenarios has very little sensitivity to the left and right selectron mass. The suppressed couplings lead to 
a reduction in the cross section. 
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V. BACKGROUND PROCESSES 



A. The Neutrino Background 



The SM radiative neutrino production is the main irreducible background for the process radiative neutrahno 
production piLl|) . The other possible backgrounds are from e+e~ — > t~^t~^, with both the r's decaying to soft 
leptons or hadrons but the contribution from this process is found to be negligible. Another large background 
is from the radiative Bhabha scattering, e+e^ — > e+e~7, where s are not detected. This radiative scattering 
is usually eliminated by imposing a cut on E^. The events are selected by imposing the condition that any 
particle other than 7 appearing in the angular range —0.95 < cos6'-y < 0.95 must have energy less than Emax, 
where Emax is detector dependent, but presumably no larger than a few GeV. This is discussed in detail in the 
literature 



The SM radiative neutrino production 



(V.l) 



has been studied extensively [43|, [SO, I69l - l71j . For this background process are produced via t-channel W 
boson exchange, and Ve,ii,T via s-channel Z boson exchange. The corresponding Feynman diagrams are shown 
inFig.Hg 
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FIG. 18: Feynman diagrams contributing to the radiative neutrino process e^e 

corresponds to the neutrinos of three flavours 



— > viyj where (iv and v) 



The background photons from this process tend to be mostly in the forward and backward directions as 
compared to the signal photons, therefore the detector acceptance cut for the photon is applied as | cos6'-y| < 
0.95. This SM background has the same photon energy distribution da/dE^ and the ^/s dependence of the 
cross section a for all the scenarios studied in this paper. Fig. [12] shows that the photon energy distribution 
from the radiative neutrino production peaks at Ej = {s — m|o)/(2Y^) « 244 GeV because of the radiative 
Z'^ production(-y/s > mz«)- The result is presented here with and without the higher order QED radiative 
effects. The inclusion of radiative effects make the peak due to the radiative return of Z'^ slightly broad. By 
imposing an upper cut on the photon energy x"^'^^ = E^'^'^/Ehcam = l~™^o/£'? , see Eq. (jIV.ll) . the photon 

background from radiative neutrino production can be reduced, by elimination of the on-shell Z'^ contribution 
to the background cross section. From the Feynman diagrams for the background process. Fig. [THl it can 
be seen that this process has a strong polarization dependence on the initial beams due to the exchange of 
W bosons which couples only to left handed electron and right handed positron. Therefore, with a suitable 
choice of beam polarization along with various kinematical cuts discussed earlier, the contribution from the 
background process can be significantly reduced. 

In Fig. [2D] we show the ^/s dependence of the total radiative neutrino cross section, with and without the 
inclusion of higher order QED radiative effects. Without the upper cut on the photon energy x™'*'', the back- 
ground cross section from radiative neutrino production e"'"e~ — vvj is much larger than the corresponding 
cross section with the cut, near the Z° production threshold. The main purpose of the cut is to move away 
from the Z'^ peak. When we impose the cut, the signal cross section from radiative neutralino production is 
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approximately three orders of magnitude smaller than the background in the case of MSSM EWSB and the 
various GUT models. 
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FIG. 19: Plot showing the photon energy 
distribution da/dEy for the radiative 

neutrino production process e^e~ — vPj 
at ^/s — 500 GeV, with and without 
including radiative effects (R.E.). 
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FIG. 20: The total energy ^/s dependence 
of the radiative neutrino cross section with 
and without an upper cut on the photon 
energy E^, along with and without the 
consideration of radiative effects (R.E.). 



e (Pi) 
>- 



e+(p2) e t 
♦ 



I>e(fcl) 



- - 

KAAAA/7(g) 



(0 



e {pi) 
>~ 



j>e(fci) 



X2n 

— < — Kaaaa/ 



(to) 



e+(p2) 



e >r 



Z 



-AAA/< 

{vii) 



7(9) 

Vl{k2) 



e (Pi 



<VW\A/7(9) 



5+^2) 



-i>e{ki) 
-K{k2) 



e (Pi) 
>- 



-fW\AA/7(g) 



e+(p2) ^2^1 
< — 



i>e(fci) 
V*e{k2) 



AAA/< 



e+(p2) 



e u 



(via) 



n{ki) 
7(9) 



e (pi) 



>- 

+ (P2) 


v+ 


<■ 


{iii 


-(Pl) 






+ (P2) 
«■ 


X2n 
X2n 



Veiki) 



WW\/7(9) 



Ueikl) 



(vi) 



FIG. 21: Feynman diagrams contributing to the radiative sneutrino production process e'^e — vv*^, with 
the last two diagrams {vii and viii) corresponding to all the leptonic sneutrino 



B. The Supersymmetric Background 

The radiative neutralino production (jlll.ip has also a supersymmetric background coming from the sneutrino 
production process (T^ 
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e+ + ^ + +7, ^ = e,/^,r. (V.2) 

This is in addition to the background from the SM process (IV.ll) . The lowest order Feynman diagrams con- 
tributing to the process (jV.2p are shown in Fig.[5T] This background process receives t-channel contributions via 
virtual charginos for I'ei^e production, as well as s-channel contributions from Z boson exchange for t'e./i.r^'e /j r 
production. In Fig. [21] we show the photon energy distribution da/dEj for radiative sneutrino production 
e+e" — j>z>*7 at y/s = 500 GeV for the MSSM EWSB and grand unified scenarios. The total cross section for 
the radiative sneutrino production is shown in Fig. [23l Since i/eV* production receives t-channel contributions 
via virtual charginos, the production cross section as well as the photon energy distribution depends on the 
chargino mixing matrix U for the different scenarios. The cuts applied to this process are the same as discussed 
before. 

As can be seen from Fig. 1231 radiative sneutrino production (IV. 21) can be a major supersymmetric background 
to neutralino production (| III.ip if sneutrinos decay invisibly, e.g. via i> — > x^v. This scenario has been called 
"virtual LSP" scenario [4J|. However, if kinematically allowed, other visible decay channels like v — x?^^ 
reduce the background rate from radiative sneutrino production. In the scenarios with a bino type neutralino, 
as in [SU{5)]24, the branching ratio BK{i/e — >■ Xi^^e) = 100%, and, therefore, this process serves as a dominant 
background. However, for the other cases where higgsino is the dominant component we have the branching 
ratios as given in Table IXIVI The second lightest neutralino being heavier than the sneutrino i>e in MSSM 
EWSB scenario, the branching ratio is kinematically not accessible. 



Branching Ratios 


MSSM EWSB 


5(7(5)75 


5(7(5)200 


SO(10)2io 


SO(10)770 


50(10)770' 


BR(u, xVe) 


78.4% 


8.1% 


21.2% 


18% 


24.2% 


44.4% 


BR(i>e ^ xW) 




1.8% 


4.54% 


0.8% 


1.2% 


6% 




21.6% 


90.1% 


74.3% 


81% 


74.8% 


49.6% 



TABLE XIV: Branching ratios of the sneutrino for different models with a higgsino type lightest neutralino 



Furthermore, neutralino production e+e~ X1X2 followed by subsequent radiative neutralino decay [73| 
X2 Xi7 is also a potential background. However, significant branching ratios BR(x2 — >■ Xil) > 10% are only 
obtained for small values of tan/3 < 5 or Mi ^ M2 |43 . ItI . [75l |. Thus, we neglect this background, detailed 
discussions of which can be found in Refs. [tJ 




FIG. 22: Plot showing the photon energy 
distribution da/dEy for the radiative 



sneutrino production process e^e 



at a/s — 500 GeV, with the inclusion of 
radiative effects. 



FIG. 23: The total energy ^/s dependence 
of the radiative sneutrino cross section 
e^e~ — ?> vv*^ with an upper cut on the 
photon energy E-^ and the inclusion of 
radiative effects. 
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C. Theoretical Significance 

Finally we address the issue of whether the photons coming from the signal process can be measured over 
the photons coming from the background SM process. The excess of signal photons Ng — crC over the SM 
background photons A'b = (JbC for a given integrated luminosity C can be expressed in terms of the theoretical 
significance [66[ 

^ =-^^. (V.3) 

A theoretical significance of S* = 1 means that the signal can be measured at a 68 % confidence level, whereas 
one needs a significance of 5 for the detection of the signal. Both the signal and the background process 
depends significantly on the beam energy for y/s = 500 GeV and £ = 5 x 10^ fb~^. In Fig. [M] we show the 
jjL dependence of the theoretical significance S for the different models considered here. When the lightest 
neutralino is dominantly a bino type state, then for /i in the range /i 6 [120,160] GeV, the significance can be a 
maximum of about 2, for the given luminosity. Then the signal for the radiative neutralino production would 
be seen at ILC as shown in Fig. [M] But for the scenarios with the higgsino as the dominant component of the 
lightest neutralino, it will be difficult to observe the signal. 
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FIG. 24: Plot showing the theoretical 
significance S for the radiative neutralino 
production as a function of /i for different 

models considered in this paper with 
= 500 GeV. 
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FIG. 25: The theoretical significance S for 

the radiative neutralino production as a 
function of the gaugino mass parameter M2 
for the different models with 
J~s = 500 GeV. 



We have also studied the variation of theoretical significance 5* as a function of the gaugino mass parameter 
M2 as well. The M2 dependence of S for all the models considered in this paper is shown in Fig. [55] in the 
interval M2 6 [200,1000] GeV. The behaviour is similar to fi with the models having a relatively higher value 
of S for lower value of M2. The values of S given here can be considered as a good guideline, since we do not 
include a detector simulation here. Besides the theoretical significance, one must also consider the signal to 
background ratio Ns/Nb in order to judge the reliability of the analysis. Overall the process under study for 
the different scenarios considered here, specially the one with a higgsino type neutralino, will most probably 
not be useful for extending the SUSY parameter space reach of ILC. However, it may be possible to reduce the 
SM background if the electron and positron beams are polarized, with right handed electrons and left handed 
positrons. 

The theoretical significance for most of the scenarios thus considered is too small, making it difficult to test 
them in the future linear colliders through the radiative neutralino production. But these scenarios have a 
distinctive feature which make it possible to test them, by means of a hard photon tag as employed for the 
process ( II.1|) . 

For some of the scenarios considered here, namely [SU{5)]i5, [SU{5)]2oo, ['S'0(10)]2io and [S'O(10)]770', due 
to large values of Mi. 2 and a low value of fi, the states Xii X2 and Xi^ are nearly degenerate, with mass around 
H, and all of them having dominant higgsino component. Since they are all closely degenerate in mass, with 
m^o — m^o around 15 GeV, therefore the processes (a) e+e^ — ?> XiX27 and (b) e+e~ X2X27 ^i^l also act 
as supersymmetric background to the radiative neutralino production. The radiative sneutrino production in 
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these scenarios with a higgsino hke LSP has visible decay channels, and can, therefore, be easily discriminated. 
Since the ZX1X2 coupling is much larger as compared to ZxiXi coupling, therefore a detailed study of the 
signatures with a hard photon and large missing energy in the final state will include processes (a), (b) and 
the radiative neutralino production. If investigated through the above channel the signal for these scenarios 
will be comparable to the SM irreducible background, in contrast to the case when only the process ( II. ip is 
considered. We show in Figs. [M] and [57] the photon energy distribution and the cross section for these models 
for different cases as follows: 

• Case 1: 

Selecting events with 7+ missing energy in the final state, which includes the process (a), (b) and the 
radiative neutralino production. 

• Case 2: 

The contribution from SM irreducible background, the radiative neutrino production. 

• For comparison we have also shown the case of radiative neutralino production XiXi7 in these Figs. 

The radiative effects are included in all these calculations. It can be seen from these Figs, that the selection 
of events with 7+ missing energy in the final state gives almost the same distribution and cross section for all 
the models considered. The models can not be discriminated but their signature is strong compared to the 
previous analyses in subsection. IIV Al The theoretical significance in this mode at y/s — 500 GeV, increases to 
about 22, compared to 0.001 for [SU{5)]j5, for the parameter values given in Table ITVl Similar result holds for 
the other scenarios. 
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VI. SUMMARY AND CONCLUSIONS 



We have carried out a detailed analysis of the radiative neutralino production e+e^ — > XiXi7 in various 
GUT models for the International Linear Collider energies and compared it with the corresponding results 
in the MSSM with universal gaugino mass parameters. In these models the boundary conditions on the soft 
gaugino mass parameters can be nonuniversal and hence different from those of MSSM with universal boundary 
conditions. This process has a signature of a high energy photon and missing energy. We have obtained a 
typical set of parameter values by excluding certain regions of the parameter space which follow from theoretical 
and experimental constraints, as discussed in Appendix. [BJ 

Using this parameter set, we have studied in detail the signal cross section for the ILC energies with un- 
polarized e"*" and e~ beams. For comparison, we have used the MSSM EWSB scenario as a benchmark. The 
contributions from the SM background e^e~ — > vv^, as well as the supersymmetric process e^e^ — !■ vv*^ 
acting as a background to the radiative neutralino production are also considered. All these processes have a 
signature of a highly energetic photon with missing energy. The photon energy distribution da/dE^, and the 
total cross section as a function of the total centre of mass energy have been calculated for MSSM EWSB, as 
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well as for the different scenarios considered here, at ^/s = 500 GeV using the programme CalcHEP. Since ISR 
and beamstrahlung are a part of the future linear colliders, due to the planned high luminosity, for a realistic 
prediction we have also included these radiative effects. The dependence of the cross section for radiative neu- 
tralino production on the SU{2)l gaugino mass parameter M2 and the Higgs(ino) mass parameter /i, as well 
as its dependence on the selectron (e^, cl) masses has also been studied and compared with the corresponding 
results in MSSM EWSB. We have considered scenarios based on SU{5) and SO{10) grand unified theories 
with nonuniversal gaugino mass at the grand unified scale. In the case of SU{5) and 5*0(10) models we evolve 
the parameters to the electroweak scale and then use these to evaluate the radiative neutralino cross section. 
All the results mainly depend on the composition of the lightest neutralino in the different models considered 
here. The models with bino as a dominant component of the lightest neutralino behave differently from the 
models with higgsino as the dominant component of the lightest neutralino. The composition of the lightest 
neutralino depends on the ratio of the soft gaugino mass parameters at the electroweak scale and fj,. The 
values for these parameters are chosen satisfying all the experimental constraints, along with the requirement 
that the lightest neutralino is the LSP. The supersymmetric background coming from the radiative sneutrino 
production is also calculated which depends on the chargino mixing matrix. The dependence of cross section 
for these two distinct scenarios is valid for different ranges of the parameters M2 and fi. Finally, in order to 
understand whether an excess of signal photons, Ns, can be measured over the background photons, N-q, from 
radiative neutrino production, we have analysed the theoretical statistical significance S = Ns/^/Ns + Nb, and 
studied its dependence on the independent parameters M2 and /z, that enter the neutralino mass matrix. It 
is seen that the signal for some scenarios is too weak to be seen at ILC. Therefore for some of these scenarios 
which also predict the second lightest neutralino to be degenerate with the LSP, the observance of a signature 
of a highly energetic photon with missing energy taking into account the second lightest neutralino production 
along with LSP production is also considered. We have also noted that initial beam polarization may reduce 
the background, and it may be interesting to study whether the signal for radiative neutralino production can 
be enhanced by using polarized beams. This question in the context of the GUT models will be studied in a 
separate paper flT\ . 



VII. ACKNOWLEDGEMENTS 



The authors would like to thank B. Ananthanarayan for many useful discussions. P. N. P. would like 
to thank the Centre for High Energy Physics, Indian Institute of Science, Bangalore for hospitality while 
this work was initiated. The work of P. N. P. is supported by the J. C. Bose National Fellowship of the 
Department of Science and Technology, and by the Council of Scientific and Industrial Research, India under 
the project No. (03)(1220)/12/EMR-II. P. N. P would like to thank the Inter-University Centre for Astronomy 
and Astrophysics, Pune, India for hospitality where part of this work carried out. 



Appendix A: Neutralino Mass Matrix, Lagrangian and Couplings 

In this Appendix we summarize the mixing matrix for the neutralinos, and the couplings that enter our 
calculations. We recall that the neutralino mass matrix receives contribution from MSSM superpotential term 

VFmssm Aii?ii?2, (A.l) 

where Hi and H2 are the Higgs doublet chiral superfields with opposite hypercharge, and /i is the supersym- 
metric Iliggs(ino) parameter. In addition to (|A.ip . the neutralino mass matrix receives contributions from 
the interactions between gauge and matter multiplets, as well as contributions from the soft supersymmetry 
breaking masses for the SU{2)l and U{1)y gauginos. Putting together all these contributions, the neutralino 
mass matrix, in the bino, wino, higgsino basis {—iX',—iX^,'ip\[_^,ip'fj^) can be written as [3. [5ll| 



■^MSSM — 



/ Ml ~mz sin dw cos /3 mz sin 9w sin (3\ 

M2 mz cosOw cos P — cos sin/3 I ( a o\ 

—mz sin 9w cos (3 mz cos 9w cos (5 ~fi I ' y ■ ) 

\ mz sin 9w sin /3 —niz cos 9w sin /3 — /i / 



where Mi and M2 are the U{1)y and the SU{2)l soft supersymmetry breaking gaugino mass parameters, 
respectively, and tan/3 — V2/V1 is the ratio of the vacuum expectation values of the neutral components of the 
two Higgs doublet fields Hi and H2, respectively. Furthermore, is the Z boson mass, and 9w is the weak 
mixing angle. We shall consider all parameters in the neutralino mass matrix to be real. In this case it can 



24 



be diagonalised by an orthogonal matrix. If one of the eigenvalues of Mmssm is negative, one can diagonalize 
this matrix using a unitary matrix iV, the neutralino mixing matrix, to get a positive semi definite diagonal 
matrix [5l| with the neutralino masses m^o {i — 1, 2, 3, 4) in order of increasing value: 

N*MmssmN^^ = diag{Tn^o, m^o, m^o, m^o) . (A.3) 

For the minimal supersymmetric standard model the interaction Lagrangian of neutralinos, electrons, selec- 
trons and Z bosons is given by [5l| 

^ = ( ^^N*,)f,PLxleR+ (iVi2 + tan0,yjVn)/ePflX?gL 

cos tiw V 2 sm 9w 



e 



4 sin 9iY cos 9w 



+eZ^/e7'' [ . a ^ TT- (I - ^in' 0w] Pl - tan OwPr] fe + h.c, (A.4) 

'- sm Ow cos &w \2 J 

with the electron, selectrons, neutralino and Z boson fields denoted by fe, eL,R, Xii and Z^, respectively, and 
Pr,l = 1 (l i 7^)- The interaction vertices following from (IA.4I) are summarized in Table IXVl 

TABLE XV: Vertices corresponding to different terms in the interaction Lagrangian (|A.4p for MSSM. Here 
we have also shown the vertices for selectron-photon and electron- photon interactions [iTj . 



Vertex Vertex Factor 

rieht selectron - electron - neutralino ""^y^ NT-, Pl 

cos -L-L ^ 

left selectron - electron - neutralino ^ „ — (V12 + tan 6wNi\)Pr 

V2 sm Oy^^- ^ ' 

neutralino - Z° - neutralino 4sin flj^cos e^.- (I^"!" " l^i*!') 7^7" 

electron - Z'^ - electron je7^ [ sin Biv^ob e-^v ~ '^''^^ ^ dwPR\ 

selectron - photon - selectron ie{p\ +P2)'' 

electron - photon - electron ie^^ 



Appendix B: Experimental Constraints 

In this Appendix we review the constraints on the SUSY particle spectrum from the different experiments 
at LHC, Tevatron and LEP. 

1. Limits on Gaugino mass parameters 

The exclusion limit on the gaugino mass parameters is set from the current experimental limits on the super- 
partner masses. Since no supersymmetric partners of the SM particles have been detected in the experiments, 
only lower limits on their masses have been obtained. In particular, the search for the lightest chargino state 
at LEP has yielded lower limits on its mass [HI. The lower limit depends on the spectrum of the model 
Assuming that toq: the soft supersymmetry breaking scalar mass, is large, the limit on the lightest chargino 
mass following from non observation of chargino pair production in e"*'e~ collisions is 

M-± > 103 GeV. (B.l) 



25 



The limit depends on the sneutrino mass. For a sneutrino mass below 200 GeV, the bound becomes weaker, 
since the production of a chargino pair becomes more rare due to the destructive interference between 7 or Z 
in the s-channel and v in the t-channel. In the models we consider, rriy is close to mo- When < 200 GeV, 
but rac, > 'm^± , the limit becomes [ssj ] 

M-± > 85 GeV. (B.2) 

For the parameters of the chargino mass matrix, the limit (|B.1|) implies an approximate lower limit (siL fssi ) 

M2, ^>100 GeV. (B.3) 

The limits shown in Eq. (jB.3[) on M2 and /i are obtained by scanning over the MSSM parameter space and are, 
therefore, expected to be model independent [12] ■ From the chargino mass limit from LEP, a neutralino mass 
limit is obtained assuming gaugino mass unification at high energy scales and amounts to 47 GeV. But due to 
the strong constraints from LHC experiments, the lower limit on x? in case of constrained MSSM (CMSSM) 
has risen to about 100 GeV. 



2. Exclusion limits on squarks and gluinos 

Experiments at the LHC and Tevatron being the proton-(anti) proton collider, with their higher centre of 
mass energies compared to LEP, is more sensitive for the SUSY particles carrying color charge, squarks and 
gluinos, because of QCD-mediated processes. Limits of the order of about 100 GeV, was set on the squark 
masses by LEP, but the recent hadron collider experiments have set much higher limits. 

Gluino masses below 800 GeV, are excluded by both ATLAS and CMS collaboration in the framework of 
CMSSM for all squark masses. For equal squark and gluino masses, the limit is around 1400 GeV [tsI. F79l|. 
These results are only slightly dependent on the choice of the CMSSM parameters, tan /?, Aq and /i. Similar 
analyses has been carried out in the context of simplified models, where upper limits on gluino pair production 
are derived as a function of the gluino and the lightest neutralino mass. For massless LSP, rrig < 900 GeV is 
excluded. For a heavy LSP above 300 GeV, no general limit on gluino can be set. 

Limits on the first two generation squark masses are set by the LHC experiments, with lowers limits of 
around 1300 GeV for all gluino masses in the framework of CMSSM. The analyses for squarks are 

similar to the gluinos. For massless neutralino in the framework of simplified models, squark masses below 750 
GeV are excluded, whereas increasing the mass of LSP above 200 GeV leads to a degradation of the limits. 

The limits on the third generation ti mass from LEP was around 96 GeV, in the charm plus neutralino 
final state 52]. LHC along with Tevatron have performed the analyses for third generation squarks in different 
scenarios, leading to different final states [soj . Similar analyses has been carried out for the sbottom quarks. 
Overall, for our analyses we are considering the scenario where the third generation squarks are excluded below 
a mass of about 800 GeV. 



3. Exclusion limit on slepton masses 

The strongest limit on the slepton masses come from LEP, because of its clean signature. LEP experi- 
ments [S^] have excluded sleptons for masses below 100 GeV in case of different scenarios. Similarly for the 
sneutrinos, limits are obtained from the invisible width of the Z boson along with the limits derived from the 
searches of gauginos and sleptons. They are excluded for a mass of about 94 GeV. 

Taking into account all the above constraints set by the LEP and LHC experiments, for our analyses we have 
taken nig « 1400 GeV, mass of first two generation squarks around 1300 GeV, around 1000 GeV along with 
the slepton masses around 150 GeV. 
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